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The freezing of metal nanoclusters such as gold, silver, and copper exhibits a novel structural 
evolution. The formation of the icosahedral (Ih) structure is dominant despite its energetic metasta- 
bility. The dynamical aspects of the structural transformations, which are eventually responsible 
for the kinetics, are studied by calculating free energies of gold nanoclusters. The transition barriers 
have been determined by using the umbrella sampling technique. Our calculations show that the 
formation of Ih gold nanoclusters is attributed to the lower free energy barrier from the liquid to 
the Ih phase compared to the barrier from the liquid to the face-centered-cubic crystal phase. 

PACS numbers: 61.46.+W, 36.40.Ei, 64.70.Nd 



Recently, nanosized metal clusters have been exten- 
sively studied as a fundamental element for technolog- 
ical applications, such as nanocatalysts and nanoelec- 
tronic devices |lj. Unlike bulk, at surfaces or inside of 
metal nanoclusters, atomic bonds may be cut and new 
bonds formed due to the presence of a nano-size sur- 
face or quantum effects. As a result, metal nanoclus- 
ters exhibit unique chemical and physical properties dis- 
tinct from bulk materials. For the controlled growth of 
low dimensional nanostructures, it is very important to 
understand the properties of metal nanoclusters, their 
formation from the liquid state or gas phase, and their 
chemistry. 

In particular, phase transitions of metal nanoclusters 
have attracted great interest because of novel physical be- 
havior, such as size-dependent melting point depression 
3, y| , quasimelting 0, |5| , and dynamic phase transitions 
(| in nano-size regions. For instance, small metal clusters 
have lower melting points compared to the bulk melting 
point [2|, |3|. Furthermore, the structure of metal nan- 
oclusters may fluctuate below the melting point under 
external perturbations [4], |5j . The experimental evidence 
for the quasimelting of gold nanoclusters was reported 
by Iijima and Ichihashi 4] through real-time microscopic 
studies: the gold clusters change from a single crystalline 
form to a twinned crystalline form based on icosahedral 
(Ih) or decahedral (Dh) structures, and vice versa, when 
irradiated with intense electron beams. 

Another interesting aspect of phase transitions of nan- 
oclusters is the novel structural evolution of clusters 
produced from the liquid state or gas phase. Because 
of a large surface-to- volume ratio in nano-size regions, 
metal clusters exhibit various structural modifications. 
For example, for face-centered-cubic (fee) noble metals 
such as gold, silver, and copper, nanoclusters of Ih or 
Dh structures with a fivefold symmetry of noncrystallo- 
graphic atomic arrangements were dominantly formed as 



observed in high resolution electron microscopy experi- 
ments |ll la, IS| • The experimental observations indicate 
that the Ih or Dh structure should be a lower energy state 
than an fee structure, by thermodynamic principles, in 
which the clusters are assumed to adopt the energetically 
stable thermal-equilibrium structure. In contrast to the 
experimental observations, previous theoretical calcula- 
tions suggest that the minimum-energy structure for Au 
is a truncated octahedron (TOh) rather than the Ih or 
Dh structure for several hundred atom clusters and the 
Ih structure is energetically metastable even for small 
clusters (less than 100 atoms) |liJ,|TjJ]. This novel and 
important structural transition behavior during the for- 
mation of metal nanoclusters thus remains ambiguous. 

In order to understand the dynamical aspects of the 
structural transformations such as the freezing and melt- 
ing of metal nanoclusters, we have performed free en- 
ergy calculations of gold nanoclusters using the umbrella 
sampling technique (12M13J combined with molecular dy- 
namics simulations. The variation of free energies for the 
structural transformation was obtained by introducing a 
crystalline order parameter |l4J . Our free energy calcu- 
lations showed that the formation of the Ih structure is 
attributed to the lower free energy barrier from the liquid 
to the Ih phase compared to the barrier from the liquid to 
the fee crystal phase. This explains why the Ih cluster is 
frequently produced in experiments despite its energetic 
metastability. 

For the description of structural evolution of gold nan- 
oclusters during the freezing and melting, we first intro- 
duce a structural order parameter $(r Ar )that represents 
the degree of crystallinity in the system |l4| . The free en- 
ergy F of the system with a particular value of structural 
order parameter $ at temperature T can be calculated 
by F($) = const - k B TlnP(^), where fc B is the Boltz- 
mann constant and P($)d<I> is the probability to find the 
system with values of $ between $ and $ + d& |l'l Il5| . 



TABLE I: 
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0.06 0.17 


0.30 


0.57 



In the above free energy calculations, conventional sim- 
ulations [16j such as molecular dynamics (MD) simula- 
tions can be used in order to determine the probability 
P($). In this case, the MD events of the structural trans- 
formations rarely occur at the transition temperature for 
nanoclusters of over a few hundred atoms. As a result, 
the probability P($) is strongly peaked around a par- 
ticular value of <I> corresponding to a certain equilibrium 
state at temperature T, while the probability is very low 
for a high free energy state. So it is practically impossible 
to obtain the free energy as a function of $ only with the 
conventional MD simulations. This problem, however, 
can be overcome by using the non-Boltzmann sampling 
such as the umbrella sampling. In the method, a "In- 
dependent bias potential makes the sampling probability 
appreciable in $ region of high free energy [1J3, yjj UJj| • 

There are two factors that play a crucial role in the 
above biased sampling method: The first is a structural 
order parameter $ as a reaction coordinate that connects 
the initial and final states of the system, and the second 
is a ^-dependent bias potential that forces the system to 
sample a high free energy region of $ space along the 
reaction path. 

For the structural order parameter, we used the bond- 
orientational parameter Qq, introduced by Steinhardt et 
al [14j . This parameter, sensitive to the degree of ori- 
entational correlation defined by the vectors connecting 
neighbor atoms, measures the degree of crystallinity of 
the system. As shown in Table[|] in a liquid cluster, Qq is 
very small (but nonzero due to the surface shape effect), 
while in a crystalline cluster, where the bond orientations 
are correlated coherently throughout the whole region, it 
is relatively large. We here note that Qq is not a state 
function with any specific internal structure but just a 
reaction coordinate indicating the degree of crystallinity. 

For the fictitious bias potential of the umbrella sam- 
pling, we chose a harmonic function localized around a 
specific value of the order parameter <!>(= Qq). We bi- 
ased the configuration space sampling by adding this bias 
potential to the original potential energy of our system. 
By adjusting the center position of the harmonic function 
to the region of a large free energy barrier [15J Il8| , we 
could achieve relatively uniform sampling over the whole 
region of the phase transition. As a result, the transition 
between cluster phases of low and high crystallinity was 
made possible despite the substantial free energy barrier. 

We now investigate the behavior of freezing and melt- 
ing of Au nanoclusters by calculating the free energies of 
the Au nanoclusters near the melting point. In these 
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FIG. 1: Representative shapes of (a) probability distribution 
and (b) free energy (AF = F(Qe) — -Fiiquid) as a function of 
a crystalline order parameter Qe at 780 K for a 459-atom Au 
cluster. 



calculations, the probability P(&) for the free energy 
was obtained by MD simulations based on the semi- 
empirical embedded-atom method |l9j and the bias po- 
tential scheme. As an initial cluster configuration, a TOh 
cluster of 459 atoms was used. The free energy was cal- 
culated at various temperatures between 760 and 820 K 
with an interval of 10 or 20 K by using the Anderson ther- 
mostat J20( . The space of bond-order parameter Qq was 
divided into 20 windows and umbrella sampling runs were 
performed with a bias potential corresponding to each 
window. Then the self-consistent histogram method |l3q 
was used to reconstruct the overall distribution of prob- 
ability P(Qq) from the individual histogram obtained in 
different windows. The total simulation time was more 
than 4 x 10 5 Ai, where At =2.5 fs represents the time 
step for the integration of the equation of motion. To 
check the reversibility of the free energy curves, we cal- 
culated the free energies by varying the Qq parameter in 
both increasing and decreasing directions. The calcula- 
tions showed no significant hysteresis except at low tem- 
perature where the free energy changes steeply with Qq. 
This indicates that our simulation time is long enough to 
equilibrate the system with negligible hysteresis. Typical 
results of the probability and free energy calculations are 
shown in Figs, ^a) and (b), respectively. More details 
on the simulation procedures and numerical techniques 
are described elsewhere |l3, |2l| . 

Figure shows the free energy (AF = F(Qq) — Fii qui d) 
curves of the gold nanoclusters as a function of a crys- 
talline order parameter Qq. Here, the liquid state is taken 
as a reference free energy state. The free energy curves 
at various temperatures contain two or more energy min- 
ima representing the liquid and crystalline structures. In 
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FIG. 2: Free energy curves for a 459-atom Au cluster 
as a function of crystallinity (Qe) at various temperatures. 
The free energy was obtained by using the umbrella sam- 
pling method combined with atomistic simulations. The inset 
shows a close-up of the curve at 770 K indicated by the box. 



order to determine the cluster structures at the minimum 
states, the typical atomic configuration at each stage was 
captured as shown in Fig. [3] The two energy minima at 
Qe = 0.45 and 0.05 represent the crystalline state of a 
TOh fee structure and the liquid state, respectively. The 
phase transition temperature is around 780 K at which 
the free energies of solid and liquid are nearly equal. At 
high temperatures such as 820 or 800 K, the liquid state 
is a global energy-minimum, while the fee solid state is 
metastable. Our free energy calculations show that at 
820 K the initial fee TOh cluster melts with a very low 
free energy barrier (less than a few k^T). When temper- 
ature is further reduced, the relative free energy of the 
crystalline solid decreases monotonically and becomes a 
global minimum state at T < 780 K. 

Here, the free energy curves below the transition tem- 
perature are particularly interesting because they are ex- 
pected to provide physical insight into the freezing be- 
havior of the Au nanoclusters. Inspection of the curves 
gives two typical freezing paths, depending on the tem- 
perature. At small undercooling temperature of 780 K, 
the curve depicts a freezing path with a very large free 
energy barrier of 10 k^T and shows the absence of any 
significant metastable solid state of a local free-energy 
minimum. However, at more undercooling temperature 
(770 K), the free energy changes significantly and a local 
free-energy-minimum state starts to appear at around 
Qq rs 0.1 [see the inset in Fig. Ej. The typical atomic 
configuration at this state shows a fivefold symmetry on 
its surface [see Fig. Efc)]. Further investigation of this 
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FIG. 3: Representative atomic configurations of a 459-atom 
Au cluster at local free-energy-minima and transition state: 
(a) The liquid state around Q@ « 0.05, (b) the transition state 
between the liquid and the Ih state, (c) the Ih state around 
Qa ~ 0.1, and (d) the TOh fee crystalline state around Qa ~ 
0.45. (e) Calculated diffraction patterns for the structures of 
(a)-(c). 



metastable structure reveals that the cluster configura- 
tion corresponds precisely to the Ih structure only with 
some defects around the center of the cluster. At this 
temperature, the TOh fee structure is thermodynami- 
cally the most stable. Here, it is noted that the free 
energy barriers are quite different: The free energy bar- 
rier for the transition from the liquid to the metastable Ih 
cluster is very low (~ k^T), while the barrier from the liq- 
uid state to the fee crystal is as large as 7 k^T. Interest- 
ingly, the typical atomic configuration at the transition 
state between the liquid and the Ih state shows {111}- 
type facets with a fivefold symmetry on its surface [see 
Fig-Bib)], although the internal structure is rather liquid- 
like as shown in its diffraction pattern [see Fig. Efe)]. As 
an atomic-level mechanism of structural transformation, 
we propose that the formation of the ordered facets pro- 
vide a way of lowering the free energy barrier between the 
liquid and the Ih state [22| . This free energy behavior is 
more prominent with further decreasing of temperature 
(760 K) as shown in Fig. [21 Considering the fast co oling 
rate of clusters in typical experimental conditions [23j, 
this simulation result clearly shows that the cluster can 
be dynamically trapped into the metastable state of the 
Ih structure with a very low free energy barrier during 
the freezing. The very low free energy barrier for the 
formation of the Ih structure explains why the Ih cluster 
structure is dominantly formed in spite of its energetic 
metastability 24] . 

In addition to the freezing behavior, the melting of 
metal nanoclusters can be also explained by the free 
energy curves of Fig. [21 Previous MD simulations re- 
ported that gold nanoclusters initially possessing an fee 
single crystalline structure underwent a structural trans- 



formation to an Ih structure before melting 25] . This 
result suggests that below the melting temperature, the 
Ih structure should be an intermediate stable state be- 
tween an fee single crystalline structure and the liquid 
phase 25] . In contrast, our free energy curves in Fig. [21 
show that the Ih structure is not stable along the melting 
path near the melting temperature. It is noted that the 
previous MD simulations were done under the constant 
energy condition J25J . The latent heat effect of melting 
arising from the constant energy condition is expected to 
affect the structural transformation. Actually, the latent 
energy lowers the temperature by more than 80 K for 
this size of nanocluster. Such instantaneous refreezing of 
the cluster leads to the formation of the Ih local energy- 
minimum state before melting (see the free energy curve 
at 760 K of Fig. |2J). In order to confirm this argument, 
we compared MD simulations under constant energy and 
thermostat conditions. In accordance with the previous 
works, we obtained the formation of the Ih structure dur- 
ing melting under the constant energy condition, but we 
could not under the thermostat condition. 

Wang et al. |26j also studied the melting behavior of 
Ih gold nanoclusters of a few thousand atoms and found 
that the cluster surface remains ordered up to the melt- 
ing temperature, although it softens considerably with 
increased diffusion due to mobile vertex and edge atoms. 
This dynamics feature also agrees well with our obser- 
vation that the high stability of {lll}-type facets near 
the melting temperature lowers the free energy barrier 
for the transition between the liquid and the Ih cluster. 

In our simulations, different cluster sizes can affect the 
freezing and melting behavior of nanoclusters. For exam- 
ple, we examined the cluster size effect by using a 561- 
atom cluster. For this cluster, the formation of the Ih 
structure is also dominant during freezing, in accordance 
with previous works [2J, |2y, [23 . Further calculations of 
the free energies of the cluster clearly showed that the low 
free energy barrier from the liquid to the Ih phase con- 
tributes to the dominant formation of the Ih structure 
during freezing. 

In summary, the freezing and melting behavior of gold 
nanoclusters was revisited in terms of the free energy as 
a function of a crystalline order parameter. By using 
the umbrella sampling technique combined with molecu- 
lar dynamics simulations, the free energy barrier for the 
structural transition was calculated. It was found that 
the dominant formation of the Ih structure during freez- 
ing of Au nanoclusters is attributed to its low free energy 
barrier. Our free energy approach not only explains well 
the experimentally observed kinetics but also provides 
complementary information for understanding dynami- 
cal properties that are experimentally inaccessible, such 
as structural properties and formation dynamics of nano- 
size clusters. It opens a way to atomic simulations, ca- 
pable of quantitatively studying the controlled growth of 
low dimensional nanostructures. 
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